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ABSTRACT 

We have measured the H i absorption towards pulsar B0329+54 using the 
Green Bank Telescope during 18 epochs between 2002 June 30 and 2003 October 
10. Three observing epochs consisted of a continuous period of 20 hours each, 
while fifteen epochs were 1-2 hours each. We calculate the structure function 
of H i absorption variations toward the pulsar on time scales of 10 minutes to 
16 months which, using the proper motion of 95 km sec" 1 and the parallactic 
distance of 1.03 kpc measured towards B0329+54 (Brisken et al. 2002), corre- 
sponds to angular scales of 0.37 /xas to 23.8 mas and samples structures between 
0.0025 — 12.5 AU assuming H i gas half way to the pulsar and ignoring scintilla- 
tion effects. We find no evidence for any turbulent H i absorption fluctuations 
towards B0329+54, with the following upper limits on At for various absorption 
features: 0.026 at —31, —21, —18, and +4kmsec" 1 0.12 at — llkmsec -1 and 
0.055 at —1km sec -1 . 

Subject headings: pulsars: individual (B0329+54) — radio lines: ISM — ISM: 
clouds — ISM: structure — turbulence 



1. Introduction 



During the last decade, studies of H i absorption lines have revealed angular or temporal 
variations corresponding to spatial scales on the order of tens of AU. In some directions 
structure in the H i has been observed and interpreted as individual H i clouds passing across 
the line-of-sight with densities of 10 4 -10 5 cm -3 . Based on these measurements a significant 
fraction (10-15%) of the cold H i gas must be in these clouds (Frail et al. 1994). However, 
it is difficult to reconcile the high H i densities at AU-scales implied by these measurements 
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with other information about the interstellar medium (ISM), as they would be greatly out 
of pressure equilibrium and should be short lived. Heiles (1997) proposed that the observed 
H i components are formed from sheets and filaments where the large column densities are 
produced by the appropriate viewing angle. Deshpande (2000) suggests that these data have 
been misinterpreted and that a single power-law describes the distribution of cold H i in 
the ISM. It has also been suggested that the observed changes in the H i absorption profiles 
towards pulsars are the result of scintillation along with a velocity gradient in a uniform 
H i medium (Gwinn 2001). Angular power spectra from H i emission observations reveal a 
power-law distribution of H i structures on parsec scales consistent with a turbulent medium 
(e.g., Green 1993). It seems possible that the AU-scale H i fluctuations are part of the same 
turbulence that is present on parsec scales, but this has not been conclusively demonstrated. 

If AU-scale fluctuations are present and are part of the turbulent cascade seen at parsec 
scales then it may be possible to determine some of the fundamental characteristics of the 
H i gas. If the gas is purely hydrodynamical (HD) then the smallest size scale that shows 
fluctuations, called the inner scale, is larger than the "molecular" mean-free-path length 
(Tennekes & Lumley 1994, Frisch 1996). For typical situations in the ISM the H i mean- 
free-path length corresponds to scales from ~ 1 — 100 AU. Collisions of the H i with the ions 
and electrons in magneto-hydrodynamical (MHD) turbulence could create H i fluctuations 
on scales smaller than ~ AU. If we are able to measure the inner scale of the HD turbulence 
then it is possible to estimate the kinematic viscosity of the H i in the ISM. The kinematic 
viscosity is given approximately by u ~ A m f p c t h, where A m f p is the mean-free-path length 
and c th is the thermal sound speed. If the H i turbulence is MHD then sub-AU scale H i 
fluctuations would be present which would tell us that magnetic fields cannot be ignored in 
any aspect of the dynamics of interstellar gas. 

We have made H i absorption measurements towards the pulsar B0329+54 with the 
Green Bank Telescope (GBT) to measure H i fluctuations on sub-AU scales. B0329+54 has 
a parallactic distance of l.OSlo'jJkpc and a proper motion of 95^ km sec -1 (Brisken et al. 
2002). The pulsar's proper motion corresponds to scales of 0.03 AU in one day and 1 AU in 
one month for an H i cloud half way to the pulsar. 
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2. Observations 

The 100 m Green Bank telescope (GBT) of the National Radio Astronomy Observatory 1 
(NRAO) was used for this measurement. The GBT is an unblocked aperture telescope with 
a spatial resolution of 9'2 at 21 cm. The system temperature on cold sky was ~ 20 K. The 
detector was the NRAO spectral processor, an FFT spectrometer, configured to have 1024 
channels for each linear polarization with a bandwidth of 1.25 MHz, producing a spectral 
resolution of 0.26 km sec -1 per channel. The accumulation memory is 32-bit, providing good 
dynamic range. The spectral processor integration time was 14 pulse periods, approximately 
10 seconds. 

The pulsar B0329+54 is an ideal target for study of small-scale structure in cold H i. 
It is very bright so a change in opacity of 0.1 can be detected in a time less than the 
scintillation time-scale of ~ 15 minutes. Its declination is such that it can be observed by 
the GBT continuously for ~ 20 hr during which ~ 50 scintles are seen. The observations 
reported here were made in 18 separate observing sessions. Three long (~ 20 hr) sessions, 
separated by two weeks, were made to probe H i fluctuations on sub-AU scales. Fifteen short 
(~ 1 — 2 hr) observations spread over the subsequent fifteen months sample scales larger than 
an AU. 

The data were calibrated using a method similar to that described by Weisberg (1978). 
For each integration an absorption spectrum was formed by taking the difference between 
the pulsar "on" and pulsar "off" spectra. If Ty is the H i opacity for the i th spectral channel 
and the j th integration sample then 

rp /„ \ p—Hj _ Tj,j(Pom Ion) ~ Tij(p ff,f on ) 

idlPonj ' " u)i E°=r Tij(p off , u) [ ) 

where n cna n is the number of spectral channels. The symbol T denotes intensity in units 
of Kelvins where, for example, Ty(p g, £ on ) corresponds to the intensity at the frequency of 
the H i line emission when the pulsar is "off", 7]j (p fr 5 ^off ) corresponds to the intensity at 
frequencies other than the H i line emission when the pulsar is "off" , and Ty (p on ) corresponds 
to the intensity at all observed frequencies when the pulsar is "on" . The power is converted 
from detector counts to Kelvin by using a calibrated noise diode that was injected every 
pulsar cycle for 10% of the pulsar period. Each integration is then weighted and summed 
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such that 

E"int T-i /_ \ „-Ti j T ?,j (Pon.^off) 
j=l 1 ij IPonJ e 

2^7=1 T 2 

where T SJ/S . is the system temperature, (Ti(p on )) is the weighted average pulsar flux, and 
there are 7ii nt independent integration samples. 

The flux from B0329+54 varies for two reasons: 1) pulse to pulse variations due to in- 
trinsic emission variations; and 2) due to scintillation. Because B0329+54 varies in intensity 
during our 10 s integration period and can be a significant fraction of the total system tem- 
perature, some of the emission spectrum will appear in the absorption spectrum as "ghosts". 
Weisberg (1978) has shown that the measured pulsar absorption spectrum is a linear com- 
bination of Ti(p off ,£ on ) and the desired pulsar absorption spectrum. Therefore the measured 
H i emission spectrum is used to fit and remove the observed "ghosts". Simultaneously a 
polynomial model is used to remove the average pulsar flux ((Tj(p on ))), which cannot be 
independently measured, and any structure in the baseline due to instrumental effects. The 
"ghost" and polynomial fit is constrained only in spectral regions where there is no H i 
absorption. The "ghost" spectrum and the polynomial fit is then extrapolated through the 
regions with H i absorption. Since the 1.25 MHz bandwidth of these observations is a signif- 
icant fraction of the scintillation bandwidth 2 for B0329+54 (measured to be 5.1 ± 0.1 MHz 
at 1640 MHz by Minter (2001) ), the average pulsar flux contains scintillation induced struc- 
tures that vary slowly with frequency. In order to remove the baseline structure arising from 
scintillation structures and instrumental effects, a 5th order polynomial is used. 

For each integration the RMS noise was computed for the part of the spectrum without 
H i line emission for both Tj(p on ) and Tj(p off ). Since the H i line emission contributes 
significantly to the system temperature (~ 20 K) the RMS noise will be larger for frequencies 
where there is H i line emission. Thus the frequency dependent noise is given by 

tfiijCpon) = o"Tj( Pon ) (l + T™/T sys ) (3) 

where TJ HI is the H i line emission strength in the ith spectral channel. A similar expression 
is used for j( Poff )- These errors are then propagated for all subsequent calculations. 

Systematic errors are also important in the comparison of H i absorption towards pulsars 
from different epochs. The systematic errors arise from the removal of the average pulsar 



2 The scintillation bandwidth is the e _1 scale over which the pulsar's signal becomes decorrelated in 
frequency due to scintillation of the pulsar signal from a non-uniform distribution of electrons in the ISM. 
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flux, (Tj(p on )), via a polynomial fit from the results of equation 2 in order to determine the 
e~ r spectra. The actual pulsar flux will certainly not follow the value of the polynomial as 
it is extrapolated over the regions of H i absorption that were excluded from the polynomial 
fitting. This difference between the actual pulsar flux and the extrapolated polynomial 
value results in a systematic error that must be taken into account when comparing the H i 
absorption spectra between different epochs. 

In Figure 1 we show an example of how the systematic error is calculated. A frequency 
band without H i emission or absorption was observed and reduced in the same manner as the 
pulsar H i absorption data. The upper panel of Figure 1 shows the results of equation 2 and 
the polynomial fit. The polynomial fit has only been done in the shaded regions, which are 
the same channel numbers used for the polynomial fit for the H i absorption measurements. 
The bottom panel of Figure 1 shows the resulting e~ T spectrum. As can been seen in 
Figure 1, there is excess noise in the regions where the polynomial fit was extrapolated (the 
white areas in Figure 1). This systematic noise is shown by the single error bar above the 
data in the bottom panel of Figure 1. The random noise is shown by the two solid lines 
in the bottom panel of Figure 1. The systematic noise was determined by taking the mean 
deviation of the e~ r spectra away from e _r = 1 within the region where the polynomial fit 
was extrapolated (the white areas in the bottom panel of Figure 1) and then subtracting, in 
quadrature, the random noise term. The systematic noise was determined by (i) subtracting 
the polynomial obtained by omitting the H i absorption channels (shown in Figure 1) with 
a polynomial obtained by using all of the channels; and (ii) taking the rms of the difference 
between these two fitted spectra over the omitted channel ranges. 

The average systematic error found using this method amounts to ±0.005 in the e~ r 
spectra. Given that previous detections of H i fluctuations have been recently questioned 
(e.g. Johnston et al. 2003, Stanimirovic et al. 2003) we choose to be conservative in our 
error estimates. The quoted uncertainties for our measurements thus contain both random 
and systematic errors, which have been added in quadrature. We use a value of ±0.005 in 
the e~ T spectra, as determined above, for the systematic errors in all H i spectra. 

3. Results 

Figure 2 compares the absorption spectrum towards B0329±54 for two epochs, July 14, 
2002 and July 2, 2003. The solid lines above and below the difference spectrum are the la 
random noise uncertainties. The systematic error from the polynomial fitting is indicated by 
the error bar above the data. There are no significant differences between the two spectra 
shown in Figure 2. In fact, no variations have been detected over the entire 16 months of 
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observations. 

The structure function for time variations in the opacity is defined as 

D T (At,i) = ((r(t,i)-T(t + At,i)) 2 ) (4) 

where r(t,i) is the opacity at time t in the ith spectral channel. The structure function 
is actually computed by binning values of (r(t,i) — r(t + At, i)) 2 and then averaging the 
data within a given bin. The errors for the structure function are computed from the 
distribution of values within a given bin. Analysis of the data using a structure function 
has two advantages: (1) it reduces the noise for a given delay through averaging; and (2) 
it provides information on the power spectrum of H i absorption fluctuations. If the power 
spectrum of opacity variations is a power law, as is the case for turbulent fluctuations, then 
the structure function will also be a power law. If only noise is present in the data then the 
structure function will be constant at a level of 2cr^. (We refer the reader to Spangler et 
al. (1989) for a more detailed discussion of structure functions.) 

Figure 3 plots the structure function for our observations of B0329+54 for two frequency 
channels. The structure function data are consistent with there being no turbulence in the 
H i absorption down to our detection limits. Figure 4 shows our observational upper limits 
to the turbulent fluctuations in opacity derived from the structure functions in each spectral 
channel. The data cover differences between ten minutes to sixteen months corresponding 
to angular scales of 0.37 yuas to 23.8 mas and geometric linear size scales of 0.0025 — 12.5 AU, 
assuming H i gas half way to the pulsar and using the pulsar's proper motion velocity of 
95kmsec _1 and distance of 1.03 kpc. (For the remainder of this paper we calculate distance 
scales assuming a distance halfway to the pulsar of 515 pc.) Furthermore, comparing our 
spectra to those of Gordon et al. (1969), who first detected H i absorption towards B0329+54, 
shows no evidence for variations in H i absorption greater than Ae~ T ~ 0.01 (the noise level 
of the Gordon et al. measurements is Ae~ T ~ 0.01 while the noise for our measurements is 
Ae~ r ~ 0.002) with a velocity resolution of 1.7 km sec -1 . This corresponds to a linear scale 
of ~ 350 AU. 

Strong interstellar scintillation is observed towards B0329+54. The multi-path propa- 
gation that results from scintillation effectively creates a spatial smoothing of H i absorption 
structures within the observed angularly broadened size of the pulsar which might mask 
absorption variations that exist on very small scales. Semenkov et al. (2003) have limited 
the angular broadening size of B0329+54 at 1600 MHz to be < 1.8 mas corresponding to a 
scintillation smoothing scale of < 1.15 AU assuming a A 2 scaling of the angular broadening 
size. Our observations probe scales up to 12.5 AU and thus scintillation cannot account for 
a lack of H i absorption fluctuations in the data. The scintillation smoothing provides a 
lower limit to the smallest scale opacity structure that we can probe. However we only have 
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upper limits to the scattering disk size. Therefore we use the geometric linear size scale of 
0.0025 AU (see above) rather than the scattering disk size as the smallest size scale probed 
by our observations. 

Gwinn (2001) suggests that interstellar scintillation coupled with gradients in the Doppler 
velocity of H i can produce small-scale fluctuations in H i absorption spectra towards pul- 
sars. The GBT observations of B0329+54 can decouple density fluctuations from velocity 
gradients since an absorption spectrum can be measured for individual scintles. The data 
points in Figure 3 are divided into three groups. The cross symbols are data averaged over 
a single scintle, the square symbols consist of ~ 2 hr averages within a given epoch, and 
the triangle symbols are data averaged over a single observing epoch. No turbulent fluctu- 
ations are detected on any timescales from the scintillation timescale to sixteen months in 
our observations. 

Gwinn (2001) predicts that the opacity variations that would be observed are given 



er T = . where Av is the amount by which the velocity of the absorption feature 



changes and c s is the thermal sound speed. We have fit Gaussians to the peak absorption 
features for each epoch. From these fits we find no trend in the value of the line center 
and can limit any change in velocity of the absorbing gas to Av < 0.065 km sec -1 . We can 
derive upper limits to the thermal sound speed from the H i spin temperature (T s ) found 
from comparing the H i emission spectrum (T em ) with the H i absorption spectrum. The 
spin temperature is found from T s = 1 T ° n ; T and should be considered an upper limit since 
some warm H i contributes to T em but typically does not contribute to the H i absorption. 
The spin temperature upper limits are shown in Figure 5. The sound speed can then be 
estimated from c s = 0.093a/T s (K) km sec -1 . An estimate of Gwinn's a T prediction is shown 
as the dotted line in Figure 4. The derived spin temperatures are not likely to be incorrect 
by more than a factor of a few. We thus consider the predictions using Gwinn's formula for 
opacity variations shown in Figure 4 to be a reasonable upper limit for opacity variations 
induced by the combination of scintillation and velocity gradients. From Figure 4 it can be 
seen that our measured limits are smaller than the upper limits determined for Gwinn's o T 
formula. Although this suggests that Gwinn's hypothesis could be incorrect, we cannot claim 
this with any certainty since only an upper limit can be determined for Gwinn's prediction. 




4. 



Discussion 



Our GBT observations towards B0329+54 were made primarily to probe sub-AU H i 
structures. The structure functions of the change in opacity versus time are consistent with 
noise (i.e. the structure functions are constant values). The structure functions thus place 
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upper limits on the opacity variations as discussed in § 3. So we do not detect H i opacity vari- 
ations consistent with a turbulent power law distribution on scales < 12.5 AU greater than 
0.026 for the —31, —21, —18, and +4 km sec -1 absorption lines, 0.12 for the —11 km sec -1 
absorption line, and 0.055 for the —1km sec -1 absorption line. This is somewhat surprising 
since Frail et al. (1994) detected H i variations of At ~ 0.1 in all observed pulsars. However, 
Johnston et al. (2003) have made multi-epoch observations of H i absorption towards four 
southern pulsars and find no significant variations and Stanimirovic et al. (2003) reach a 
similar conclusion from re-observations of the Frail et al. pulsars to increase the number of 
temporal baselines. VLBA H i absorption measurements have been summarized by Faison 
(2002). Only two sources show significant H i fluctuations. Small-scale H i structure mea- 
sured through H i absorption is currently found in only two of the 15 sources observed in 
Johnston et al. (2003), Stanimirovic et al. (2003), Faison (2002) and this work. 

On large scales the distribution of H i is influenced by spiral density waves, supernovae, 
etc. as is evident by the observed shells and filaments in neutral hydrogen surveys. Dickey 
& Lockman (1990) have argued that no more than 10% of the total H i exists in small-scale 
structures (< lpc). That is, there are not large variations in the H i column density on 
small spatial scales and the concept of distinct H i clouds does not describe most of the 
Galactic neutral hydrogen. The distribution of H i can be analyzed by producing angular 
power spectra of H i emission over different spatial scales. The results are well fit by a 
power-law with a slope of approximately —3 towards different directions in the Galaxy (e.g., 
Crovisier & Dickey 1983; Green 1993; Dickey et al. 2001). These results can be described 
by a turbulent cascade of energies (Lazarian & Pogosyan 2000). But does this turbulence 
extend down to very small spatial scales? In other words, what is the inner scale of the 
turbulence? 

We are unable to answer this question with the current data on B0329+54 since we have 
not detected any significant turbulent H i fluctuations. 

Of particular note are the results of Shishov et al. (2003) from diffractive scintillation 
measurements of B0329+54. They find that on scales below 3 x 10 16 cm that the diffractive 
scintillation can be explained with a scattering screen comprised solely of ionized gas. 3 How- 
ever, on scales above 3 x 10 16 cm their results require that some neutral gas is also present 
within the scattering screen. Could this indicate that the inner scale for the neutral gas is 
2000 AU? If this is the inner scale then an H i density of 1 cm~ 3 at a temperature of 100 K 
would give a kinematic viscosity of 3 x 10 21 cm 2 s -1 . 



3 The value 3 x 10 16 cm was obtained by Shishov et al. (2003) using a velocity of 139 km sec 1 for the 
pulsar. Using a velocity of 95 kmscc -1 as measured by Brisken et al. (2002) results in a scale of 2 x 10 16 cm. 
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Fig. 1. — Top panel: A spectrum of (Tj(p on )) e~ n for two polarization (X and Y) towards 
0329+54 and offset in frequency from the H i line. The smooth, solid curves are the polyno- 
mial fits to the data. The grey area denotes the channels used in the polynomial fit where 
as the white area is where the polynomial was extrapolated. The same regions (by channel 
number) used here were also used for the H i line data. Bottom Panel: The resulting spec- 
trum of e~ n . The "straight lines" indicate the la level of the random noise and the "thick" 
error bars indicate the la errors in extrapolating the polynomial. 
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Fig. 2. — Top: The H i emission spectrum towards B0329+54. Middle: The absorption 
spectra from July 14, 2002 (solid line) and July 2, 2003 (dotted line). Bottom: Difference 
between the two absorption spectra. The lines above and below the spectrum are the la 
uncertainties. The uncertainty is a function of frequency since the H i emission contributes 
significantly to the total system temperature. The la uncertainty in extrapolating the 
polynomial fit to the pulsar emission spectrum across the H i line is shown by the individual 
error bar. 
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Fig. 3. — Structure function of r for the spectral channels centered on vlsr = —11.6 km/s 
(left panel), at the peak of the absorption line, and vlsr = —15.7 km/s (right panel), halfway 
down the absorption line, with 3<x error bars. The cross symbols are data averaged over a 
single scintle, the square symbols consist of ~ 2 hr averages within a given epoch, and the 
triangle symbols are data averaged over a single observing epoch. The best fit power laws 
are shown as the solid lines and are given by (—1.5 ± 0.1) + (0.00 ± 0.02) • log(At) and 
(-2.4 ± 0.1) + (-0.02 ± 0.02) • log(Ai) respectively. No turbulent H i opacity fluctuations 
are detected on any scale between 0.0025 and 12.5 AU. 
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Fig. 4. — Upper limits to At from our measurements (solid line) and from the scintillation 
model of Gwinn (2001) (dotted line). The frequency dependence of the Ar upper limits 
follows the expected noise contributions due to the change in system temperature from the 
H i emission and absorption. 
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Fig. 5. — Upper limits to the spin temperature, T s , of the H i responsible for the absorption 
towards B0329+54. 



